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Abstract. Recently there has been a huge activity in the dark matter direct detection field, 
with the report of an excess from CoGeNT and CRESST along with the annual modulated signal 
of DAMA/Libra and the strong exclusion bound from XENONIOO. We analyse these results 
within the framework of Bayesian inference and evidence. Indeed bayesian methods are well 
suited for marginalizing over experimental systematics and background. We present the results 
for spin-independent interaction on nucleus with particular attention to the low dark matter 
mass region and the compatibility between experiments. In the same vein we also investigate 
the impact of astrophysical uncertainties on the WIMP preferred parameter space within the 
class of isotropic dark matter velocity distributions. 



1. Introduction 

Dark Matter (DM) direct searches experience a fervent activity with the report of an annual 
modulation by the CoGeNT collaboration pQ along with the established DAMA/Libra (DAM A 
hereafter) signal The exclusion bound of the XENONIOO collaboration [3] (XelOO hereafter), 
obtained via a profile likelihood approach puts to test those excesses. The aim of this 
proceeding is to update the bayesian analysis of [5] with the CoGeNT data 2011 and to 
investigate the mutual agreement/disagreement between those experiments for spin- independent 
interaction. As a novelty, we present the assessment of the need of including astrophysical 
uncertainties in the direct detection (DD) analysis, with the use of the Bayesian evidence. 

2. Bayesian framework and experimental likelihoods 

The bayesian statistics is a well defined framework that allows the inclusion of experimental 
systematics and backgrounds, as well as astrophysical uncertainties common to all experiments. 
Indeed the posterior probability density function (pdf) is given by Bayes' theorem: 



where the likelihood C describes how the theoretical model, with free parameters 0, describes the 
data X. The prior probability density vr encodes the state of knowledge on Q before observing 
the data and is therefore independent of X. The posterior pdf is sampled via MCMC techniques 
using the package CosmoMC (BJ ^ and then marginalized over nuisance parameters. For DD 
purposes, the important quantity is 'Pmarg('TiDM, o"^'^), where ttidm is the DM mass and is 
the DM-nucleus cross-section. 
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We take as nuisance parameters the mean velocity of the DM at the sun position vq, the 
escape velocity from the DM halo Vesc and the DM density at the sun position pq (altogether 
named astrophysical variables). We consider two models for the DM velocity distribution: SMH 
(Standard Model Halo) - that is a simple Maxwellian halo with astrophysical variables fixed at 
their mean value - and NFW - namely we derive the velocity distribution from this DM density 
profile and marginalize over all possible values of vq, Vesc and Pq as well as profile parameters. 
For all technical details about the analysis, prior choices and likelihood definitions, we refer 
to [5j. Here we just recall that the main systematic for XelOO is the scintillation efficiency, while 
for DAMA the quenching factor on Sodium gNa and on Iodine qi are marginalized over all their 
experimental range. Regarding the new CoGeNT data, the likelihood is given by the product 
of two gaussian distributions, one for the total rate and one for the modulated signal. The 
data binning, analysis and subtraction of cosmogenic peaks follow closely [1] and [8]. A flat and 
exponential background, which does not modulate in time and which is described by 3 nuisance 
parameters, is added on top of the DM signal. The exponential background accounts for the 
effect of a bad rejection of surface events near threshold. 

In Eq. [Tjthe denominator Z{X), called Bayesian evidence, is defined as the average of the 
likelihood over the prior for a specific model Ai 

Z = J C{X\e)Tr{9)d^9 , (2) 

where D is the dimensionality of the parameter space. This quantity is of crucial importance in 
model comparison, which is defined as the ratio of posterior probabilities. A model A^o can be 
compared with Aii, which has extra parameters, through the Bayes factor B, see e.g. P [TU]: 

V{Mo\X) ^ ^o 7r(.Mo) ^ ^(.Mp) 
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where 7r(A^j) is the prior probability of each model. The Bayes factor automatically favours 
simpler models unless the data justify the complexity of more complicated alternatives, because 
of the marginalization procedure used to calculate the evidence, Eq. [2| In the following we 
will compare the model SMH with NFW, the evidence for each model being computed with 
MULTIMEST pH[T2]. 

3. Results 

Figure [l] summarizes the 2D inference in the {ttidm, c^^} plane for individual experiments 
(blue dashed curve for XelOO, gray dashed and solid black contours for CoGeNT and DAMA 
respectively) and for the combined fit (solid red lines), in the SMH and NFW model (left and 
right panel respectively) . Table [l] resumes the preferred values for the DM and astrophysical 
observables in each model (recalling that vq, Vqsc and Pq are fixed for SMH model). Comparing 
the left and right panel one notes that the closed regions become wider, while the XelOO exclusion 
bound moves slightly to the right, improuving the agreement between various search results. This 
is a volume effect due to the marginalization procedure over the astrophysical variables. In both 
models, the inclusion of uncertainties on the scintillation energy leads to compatibility between 
CoGeNT, the combined fit and XelOO exclusion bounds at 905"% C.L., while leaving a marginal 
compatibility with DAMA at 99% C.L.. It is always possible to combine experiments and find 
a common region, but actually the question is whether the preferred point for {mDM><7n^} is 
a good fit for both experiments and what are the values for the nuisance parameters that are 
selected. Notice from table [T] that the combined fit chooses values which are in line with the one 
of CoGeNT and DAMA alone, within the statistical errors. The large deviation from DAMA 
alone fit concerns the quenching factor on Sodium: as it has been shown in figure 2 of |5j the ID 




Figure 1. Left: 2D credible regions for the individual experimental bounds and regions assuming 
the SMH, combined in a single plot. For DAMA (black solid), CoGeNT (gray dashed) and the 
combined fit (red solid) we show the 90% and 99% contours. The blue dashed line represents 
the 905% confidence level (C.L.) for XelOO. Right: Same as left for the NFW model. 



posterior pdf for QNa is flat all along its prior range, while the combined fit is strongly peaked 
at the value gNa = 0.57 it 0.03. This can be understood thinking that the larger gNa the lower 
the DM mass should be to account for the DM modulated signal. The selected value for the 
combined fit therefore pushes the DAMA region towards the direction of the CoGeNT region. 

Table [2] resumes the results for model comparison in the case of individual experiments and for 
the DAMA+CoGeNT combined fit. A positive (negative) value for InB represents an increase 
(decrease) of the support in favour of the simplest model given the observed data. The strength 
of evidence is given by the empirical 'Jeffreys scale' (see table 1 in flOj), where threshold values 
are set for inconclusive, moderate or strong evidence: e.g. Ini? = 5 states that the odds against 
the most complicated models are 150:1 and corresponds to a probability of 0.993. This holds 
in the case where the probability of each model is the same, namely 7r(A^o) = tt{M.i) = 1/2. 
Looking at In B in table [2] it is clear that there is strong or moderate evidence against NFW 
model for single experiments. Therefore data at present time do not allow a determination 
of astrophysical observables. On the contrary there is a very strong evidence for NFW for 
the combined fit, namely experiments need to adjust astrophysical observable values to find a 
common agreement. Similar results hold for the whole class of spherically symmetric DM halos. 

4. Conclusion 

We have employed Bayesian inference to study the concordance between various DD experiments 
and model comparison to infer the necessity of a more complicated model than the SMH at 
present time. The outcomes of the analysis may be summarized in few main points: 

• the combined fit of DAMA and CoGeNT requires parameters closely in line with those of 
individual fits, the only exception being a large quenching factor on Sodium; 

• the CoGeNT region and the combined fit are compatible at 90^% C.L. with the XelOO 
exclusion bounds, while DAMA is only marginally allowed at 99% C.L.; 

• a combination of at least two experiments demonstrates the ability of constraining 
astrophysical variables, on the contrary of a single DD experiment. 



Table 1. ID posterior pdf modes and 90% credible intervals for the circular velocity vq, escape 
velocity Vcsc, and the local DM density pQ for the 2 DM halo model considered in this work. 



Model 



m-DM (GeV) 



(cm^) Vq (km/s) 



(km/s) Pq (GeV/cm3) 



SMH DAMA 
CoGeNT 
Combined 
XenonlOO 



12 
7.5 
6.9 



2.2 X 10 
1.9 X 10 
1.8 X 10 



-40 
-40 
-40 



230 
230 
230 
230 



544 
544 
544 
544 



220;^° 
219tl 
214^^^ 

910+43 



0.4 
0.4 
0.4 
0.4 



NFW 



DAMA 
CoGeNT 
Combined 
XenonlOO 



12 
7.4 
7. 



1.7 X 10-4° 
9.7 X 10-41 
1.5 X 10-40 



558] 
557^ 



3-1-19 
-16 
7+18 

-15 

556^14 
559 ± 18 



-1-1-0.15 
-0.09 
7-f0.15 
-0.10 

^■"^3-0.09 

7-1-0.16 
-0.08 



0.39] 
0.37^ 



0.37: 



Table 2. The odds against the NFW model for individual and combined experiments. 



DAMA CoGeNT Combined 
\nB 5.3 3.9 -32 



Remarkably the excess region falls in the same ballpark of {mDMjfn^} values that could 
explain the excess of events reported by the CRESST experiment p!3]. The CoGeNT 
collaboration has presented an ongoing analysis on the rejection of surface events near 
threshold [14], which could account for at least 50% of the counts at low energies: the total 
rate would therefore be accommodated by smaller DM cross-sections than those in figure [T} 
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